A Fermion Crystal with Quantum Coherence 
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When two-dimensional electrons are subjected to a very strong magnetic field, they are believed 
to form a triangular Wigner crystal. We demonstrate that, in the entire crystal phase, this crystal is 
very well represented by a composite-fermion-crystal wave function, revealing that it is not a simple 
Hartree-Fock crystal of electrons but an inherently quantum mechanical crystal characterized by a 
non-perturbative binding of quantized vortices to electrons, which establishes a long range quantum 
coherence in it. It is suggested that this has qualitative consequences for experiment. 

PACS numbers: 73.43.-f,71.10.Pm. 



Wigner argued that electrons form a crystal when the 
interaction energy of electrons dominates over their ki- 
netic energy |l|. O ne such situation occurs when a two- 
dimensional electron system is exposed to a strong mag- 
netic field (B). The kinetic energy is quantized into 
Landau levels under a magnetic field, and the number 
of occupied Landau levels is called the filling factor, 
v = phc/eB, p being the two-dimensional electron den- 
sity. At sufficiently strong magnetic fields, when all elec- 
trons fall into the lowest Landau level (v < 1), the kinetic 
energy is no longer relevant, and the nature of the state 
is determined solely by the Coulomb interaction. One 
might expect a crystal for any v < 1, but for a range of 
filling factors the system condenses into a quantum Hall 
liquid, characterized by dissipationless transport and pre- 
cisely quantized plateaus of Hall resistance |2j • There are 
strong indications, however, that a crystal occurs at suf- 
ficiently low v, and its properties have been studied in 
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Our understanding of the physics of the FQHE liq- 
uid is secure in part because the microscopic wave func- 
tions |24i J2R| incorporating this physics have been shown 
to be extremely accurate through comparisons against 
exact wave functions known for systems containing a fi- 
nite number of electrons. The liquid wave functions be- 
come progressively worse with decreasing v, consistent 
with the expectation that a crystal becomes energeti- 
cally favorable at some point. Indeed, a Hartree-Fock 
wave function describing an "electron crystal" (EC) pro- 
vides a better approximation for the ground state be- 
low certain filling factor. However, the electron crys- 
tal is not particularly good quantitatively. The aim of 
this paper is to show that the low v phase is strikingly 
accurately described as a quantum crystal of composite 
fermions, where composite fermions are bound states of 
electrons and quantized vortices. The composite- fermion 
(CF) crystal wave function offers a quantitative descrip- 
tion of the crystal state at a similar level of accuracy 
as the CF liquid wave functions do for the FQHE liquid. 
One might have expected residual CF correlations to sur- 



vive in the crystal phase in the proximity of the phase 
boundary separating the CF liquid and the crystal, but 
we find that the CF crystal is realized even deep inside 
the crystal phase, down to arbitrarily low fillings. We 
estimate the temperature below which thermal fluctua- 
tions should not destroy the quantum effects due to the 
formation of composite fermions in the crystal. 

The gaussianly localized wave packet for an electron in 
the lowest Landau level localized at R = (X. Y) is given 
by |H 
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where the magnetic length, Iq = ^Jfic/eB, has been taken 
as the unit of length. The wave function for the elec- 
tron crystal (EC) is constructed by placing electrons on 
a triangular lattice Rj , the lowest energy solution for the 
classical problem, and then antisymmetrizing the prod- 
uct: El 
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where the sum is over all permutations P and ep is 
+1 for even permutations and —1 for odd permuta- 
tions. With the lattice constant a = (A'k/\/Zv) 1 / 2 Iq, the 
overlap integral between nearest neighbor electron wave 
functions [Lj, exp(— a 2 /2^) = exp(— 3.627/ 1/), decays 
rapidly with decreasing v. We will work with the sym- 
metric gauge, A = (B/2)(—y,x,0), for which the total 
angular momentum L is a good quantum number. Be- 
cause the wave function ty EC is not an eigenstate of an- 
gular momentum, we follow the method of Yannouleas 
and Landman 26] to project it onto a definite L, denot- 
ing the resulting wave function ty^. Such projection 
amounts to creating a rotating crystal, implying that the 
crystalline structure is not apparent in the density but in 
the pair correlation function. The explicit expression for 
*f c is given in Ref. H(| 

Following the standard procedure of the composite 
fermion (CF) theory |25l.l27|. we construct the following 
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FIG. 1: The correlation energy of the optimal CF crystal, 
i.e., the % deviation of its Coulomb energy from the Coulomb 
energy of the uncorrelated electron crystal, for N = 6 par- 
ticles. The superscript 2p on 2p CFC indicates the vortex 
quantum number of composite fermions. The energy of the 
electron crystal for L > 400 is taken from Yannouleas and 
Landman |2ql . The deviation of the exact energy from the 
electron crystal energy is also shown for L < 145; for larger 
angular momenta, where the exact energy is not available, 
we show an accurate approximation, VqJ, (explained in the 
text), as an independent reference. For 2p > 6, the number 
of vortices carried by composite fermions is shown in brackets 
near the diamond. The energy difference per particle between 
the electron and the CF crystals is given in the inset, quoted 
in units of e 2 /eZo, where lo is the magnetic length and e the 
dielectric constant of the host semiconductor. 



wave function: 
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It is interpreted as a CF crystal (CFC), because the Jas- 
trow factor n ? <fc( z j — Zk) 2p binds 2p quantized vortices 
to each electron in ^ EC to convert it into a composite 
fcrmion. The composite fermions of different flavors are 
denoted by 2p CF, and their crystals by 2p CFC. A vari- 
ational wave function analogous to that in Eq. © was 
considered earlier by Yi and Fertig 20] and was shown to 
have lower energy than the Lam-Girvin 15| wave func- 
tion in the filling factor range 0.1 < v < 0.2. Narevich, 
Murthy, and Fertig |2lJ used a Hamiltonian formulation 
of composite fermions to estimate gaps and shear mod- 
ulus on either side of the v — 1/5 quantum Hall state. 
However, it is not known how close ty^ CFC is to the 
actual ground state and what is the filling factor range 
of its applicability. It is noted that Jastrow-type wave 
functions were earlier used to estimate the effect of zero 
point motion in the 4 He crystal 28]. 



Comparisons with exact solutions have served as a rig- 
orous and unbiased means for establishing the physics of 
various strongly correlated states in the lowest Landau 
level. We now subject ^ PCFC to such tests. Exact wave 
functions can be obtained (using the Lanczos method) 
for up to N = 7 particles in the low-^ region of interest. 
We will present below detailed results for N — 6; the 
study of N = 5 and N = 7 particles is consistent with 
our conclusions below. The filling factor of the finite sys- 
tem will be defined by the expression v = N(N — 1)/2L, 
which gives the correct value of v in the thermodynamic 
limit. For N = 6, the lowest energy classical configura- 
tion has one particle at the center, with the remaining 
five forming a ring around it 29] . The wave functions in 
Eq. IJ2J) for 2p ^ have rather complicated correlations 
built into them, but the interaction energy per particle, 
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can be evaluated by the Metropolis Monte Carlo method 
at least for many large values of 2p (the computation 
time increases rapidly as 2p is reduced). The total en- 
ergy also has contributions from electron-background and 
background-background interactions, but these terms are 
the same for different crystal wave functions for a given 
L, so are not relevant for comparisons. 

Fig- El shows the correlation energy of the optimal CF 
crystal, defined as the deviation of its energy from that of 
the uncorrelated electron crystal. (The computationally 
accessible range of 2p allows us to determine the mini- 
mum CFC energy. The only exceptions are the largest 
two values of L, where we show the energy at the smallest 
2p studied; the minimum CF energy may be still lower 
here.) The exact correlation energy is also shown for 
L < 145. For L > 145, the dimension of the Fock 
space (D) is too large for an exact treatment. As an 
independent reference point, we obtain an accurate ap- 
proximation to the exact energy by the method of "CF 
diagonalization," wherein the Coulomb Hamiltonian is 
diagonalized in a correlated CF basis, the dimension of 
which is much smaller than the dimension of the full basis 
needed for the exact state; gradually increasing the ba- 
sis size gives an increasingly better approximation [3fJ. 
Fig. ^ quotes Vqp (using the notation in Ref. l30j) , ob- 
tained with a correlated CF basis of dimension 150. V, 



(2) 
CF 



has been shown to be very precise 30] : for the six particle 
system it is within 0.02% of the exact energy for L < 145 
and we expect similar level of accuracy for higher L as 
well. 

The minimum energy for all L is obtained at a non-zero 
value of 2p, which establishes that the CF crystal pro- 
vides a better variational state than the electron crystal. 
Most significantly, the CF crystal is essentially the exact 
solution for v < 1/7 (L > 105). For 100 < L < 145, 
the energy of the optimal CF crystal is approximately 



TABLE I: The last three columns give the overlaps of 
CF crystal (CFC), electron crystal (EC), and Laughlin's 
wave function with the exact ground state wave function 
at several filling factors v. The overlap is defined as 

u^trial\^exact\l2 , ,^triall^,trial\ /^,exactl^,exact\ rj-^e second 

column gives D, the dimension of the basis space for N — 6 
electrons, and L is the total angular momentum of the state. 



v (L) 


D 


CFC 


EC 


Laughlin 


1/5 (75) 


19858 


0.891 


0.645 


0.701 


1/7 (105) 


117788 


0.994 


0.723 


0.504 


1/9 (135) 


436140 


0.988 


0.740 


0.442 



within 0.02% of the exact energy. Tables HI and ITT1 show 
how Laughlin's liquid wave function, *f> EC , and ty CFC 
compare with the exact wave function for v = 1/5, 1/7 
and 1/9. As indicated earlier, the liquid wave function 
worsens and ^> EC improves with decreasing v, but nei- 
ther is very good at small fillings. In contrast, *B CFC is 
surprisingly close to the exact state. Its overlap with the 
exact wave function is ~99% for v = 1/7 and 1/9, while 
its energy deviates from the exact energy by 0.016% and 
0.006%, respectively. To appreciate the significance of 
these comparisons, we note: (i) The exact state is a lin- 
ear combination of a large number of Slater determinant 
basis functions (see Table UJ, involving D — 1 parame- 
ters, and yet, a single CFC wave function captures its 
physics almost exactly, (ii) The CFC wave function for 
v < 1/7 is more accurate than accepted wave functions 
for the FQHE liquid states. For example, for the sys- 
tem of six particles at v = 1/3 in the disk geometry, the 
energy of Laughlin's wave function J24J is off by 0.15% 
and its overlap with the exact state is 0.964, in spite of 
the fact that the dimension of the Fock space is much 
smaller (D — 1206). For larger L, the energy of the CFC 

(2) 

is lower than V^p , with the possible exception of the last 
two points, where we may not have the optimal CFC. We 
take these comparisons as a definitive verification of the 
formation of a CF crystal at low v. 

Because every particle sees quantized vortices on ev- 
ery other particle, the formation of composite fermions 
implies a long range quantum coherence in the crystal 
phase. To get a feel for how the binding of vortices to 
electrons affects the inter-particle correlations, we show 
in Fig.|2the pair correlation function g[x) for several can- 
didate wave functions as well as the exact ground state 
for v = 1/7; g(x) is the probability of finding a pair of 
particles at an arc distance ions circle of radius R. 
(R is chosen to match the distance of a particle in the 
parent classical crystal from the center of the disk.) The 
result shows that the crystalline correlations are slightly 
weakened by the formation of composite fermions. It is 
perhaps counter-intuitive that such an effect should lead 
to a lower energy even at very low fillings. 



TABLE II: Interaction energies per particle for the exact 
ground state, the CF crystal (CFC), the electron crystal (EC), 
and Laughlin's wave function for six particles at several filling 
factors. The uncertainty in the last digit from Monte Carlo 
sampling is given in parentheses. 



v (L) 


exact 


CFC 


EC 


Laughlin 


1/5 (75) 


2.2019 


2.2042(5) 


2.2196 


2.2093(2) 


1/7 (105) 


1.8533 


1.8536(2) 


1.8622 


1.8617(2) 


1/9 (135) 


1.6305 


1.6306(1) 


1.6361 


1.6388(1) 



Of interest is the nature of the thermodynamic state, 
obtained in the limit N — > oo at a fixed filling factor. 
Finite size studies do not necessarily provide a reliable 
insight into the thermodynamic state. For example, at 
v = 1/5 the CFC gives the best description for N = 6 
but the thermodynamic state here is known to be a liq- 
uid pj. |_lJ] . However, when an extremely precise and un- 
ambiguous description of the finite N state is obtained, 
as is the case at v < 1/7, we consider that to be a strong 
indication for the nature of the state in the thermody- 
namic limit. In any case, even though our finite N study 
cannot give the precise v value where a transition from 
liquid to crystal takes place, it does make a compelling 
case that whenever the thermodynamic state is a crystal, 
it is a crystal of composite fermions, even in regions of 
the phase diagram far from the CF liquid. 

The quantum character of the crystal is not fragile, 
and ought to be observable at presently attainable tem- 
peratures, even at very small v. The energy difference 
per particle, V CFC — V , shown in the inset of Fig. ^ 
gives a crude estimate for the temperature below which 
the quantum nature of the crystal should be robust to 
thermal fluctuations. The relevant temperatures appear 
to be well within the present experimental reach - for 
example, for parameters appropriate for GaAs, the quan- 
tum crystal regime is estimated to be below sa 25 mK (at 
B = 25T) even at v = 1/33. From the N dependence, we 
have estimated that the energy difference shown in the 
inset underestimates the thermodynamic energy differ- 
ence by approximately a factor of two. It is interesting 
to note that even as the energy difference between the 
CF and the electron crystals decreases as v — * 0, 2p con- 
tinues to rise. Thus, CF flavors of up to very high 2p 
are predicted to occur in the crystal state. In the liq- 
uid phase, 2 CFs and 4 CFs have definitely been observed, 
and there is evidence also for 6 CFs and 8 CFs at relatively 
high temperatures, |l(j 

Quantum crystals of neutral bosons have been studied 
previously, where quantum mechanical exchange due to a 
significant overlap between the neighboring atoms is re- 
sponsible for the quantum behavior, including a possible 
"supersolid" phase in which Bose-Einstein condensation 
can coexist with crystallinity |3jJ, |32( . In the lowest Lan- 







activation energy in the crystalline state, compared to 
theoretical predictions based on an electron crystal, as 
well as its non-monotonic filling factor dependence has 
also been rationalized in terms of a CF crystal J2(J, |2lJ . 

In summary, we have shown that the two-dimensional 
electron system in high magnetic fields provides an ex- 
ample of an intrinsically quantum mechanical crystal of 
fermions with non-trivial correlations. Partial support of 
this research by the National Science Foundation under 
grant no. DMR-0240458 is gratefully acknowledged. 



FIG. 2: The pair correlation functions for the CF crys- 
tal (solid circles), the electron crystal (empty squares), and 
Laughlin's wave function (empty triangles) on a circle of ra- 
dius R = 6.445/o for six particles at v — 1/7. The solid line 
shows the exact pair correlation function. 



dau level, the overlap between (uncorrelated) electron 
wave packets at neighboring sites is negligible in the fill- 
ing factor region of interest (the overlap integral is 10~ 15 
for v = 1/9). The quantum nature of the CF crystal 
owes its origin to the long range Coulomb interaction. 

Given that the CF liquid behaves qualitatively differ- 
ently from an electron liquid, one may ask in what ways 
the properties of the CF crystal are distinct from those 
of an electron crystal. We mention here a few examples 
where the CFC can provide natural explanations for cer- 
tain experimental facts, although further work will be 
needed to make the connection with experiment more 
direct and to clarify other possible implications. The 
issue is obviously relevant to experiments that exhibit 
transitions between the liquid and crystal phases. Re- 
entrant transitions between the FQHE liquid and an in- 
sulating state, thought to be a pinned crystal, have been 
observed [3, |4| in going from v = 2/5 to v — 1/5. In 
the filling factor range 1/5 > v > 1/9, the low tem- 
perature insulating state melts into a CF liquid upon 
raising temperature jl(l | , as indicated by the appearance 
of FQHE like structure. These observations become less 
baffling knowing that the crystal is itself made of com- 
posite fermions rather than electrons, thus requiring a 
less drastic reorganization of the state at the transition. 
Another result, perhaps puzzling for an electron crys- 
tal, is that the Hall resistance of the pinned crystal is 
close to the value it would have for a liquid [23| ■ If the 
current is carried by composite fermions instead, then 
the Hall voltage induced by the accompanying vortex 
current (the vortices effectively behave as magnetic flux 
quanta |25j,|24(), through an effective Faraday effect, is 
rou ghly consistent with the observation. (Zheng and Fer- 
tig |l9j considered a similar mechanism for transport by 
correlated interstitial defects.) The unexpectedly small 
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